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ABSTRACT 

We use HST/WFC3 imaging from the CANDELS multicycle treasury survey, in conjunction with the Sloan 
Digital Sky Survey, to explore the evolution of galactic structure for galaxies with stellar masses > 3 x 1O 1O M 0 
from z- 2.2 to the present epoch, a time span of 10 Gyr. We explore the relationship between rest-frame 
optical color, stellar mass, star formation activity and the structural parameters of galaxies as determined from 
parametric fits to the surface brightness profiles of galaxies. We confirm the dramatic evolution from z= 2.2 to 
thepresentday in the number density of non-star-forming galaxies above 3 x 1O ,O M 0 reported by other authors. 

We find that the vast majority of these quiescent systems have concentrated light profiles, as parameterized by 
the SSrsic index, and the population of concentrated galaxies grows similarly rapidly. We examine the joint 
distribution of star formation activity, S6rsic index, stellar mass, mass divided by radius (a proxy for velocity 
dispersion), and stellar surface density. Quiescence correlates poorly with stellar mass at all z < 2.2 (given 
the < 0.2 dex scatter between halo mass and stellar mass at z ~ 0 inferred by More et al. 2009; MNRAS, 

392, 801, this argues against halo mass being the only factor determining quiescence). Quiescence correlates 
better with S6rsic index, 4 velocity dispersion’ and stellar surface density, where S6rsic index correlates the best 
(increasingly so at lower redshift). Yet, there is significant scatter between quiescence and galaxy structure: 
while the vast majority of quiescent galaxies have prominent bulges, many of them have significant disks, and 
a number of bulge-dominated galaxies have significant star formation. Noting the rarity of quiescent galaxies 
without prominent bulges, we argue that a prominent bulge (and, perhaps by association, a supermassive black 
hole) is a necessary but not sufficient condition for quenching star formation on galactic scales over the last 
10 Gyr; such a result is qualitatively consistent with the expectations of the AGN feedback paradigm. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: structure — galaxies: evolution — galax- 
ies: general 


1. INTRODUCTION 

The last decade of study has brought into sharper focus 

than ever the bimodality of the star formation histories of 

galaxies. There is a relatively tight distribution of star for- 

mation rates (SFRs) at a given mass for star forming galax- 

ies 0.3 dex scatter, with a few outliers to high SFR; 

e.g., Brinchmann et al. 2004, Salim etal. 2007), persisting 

out to z > 2 (Noeske et al. 2007; Wuyts et al. 2011a). The 

red sequence has SFRs substantially below those expected 

for star forming galaxies (but often with some star forma- 

tion; see, e.g., Yietal. 2005). We will call these quies- 

cent galaxies in what follows. The relative prominence of 

the two populations is a function of stellar mass, surface 


density, inferred velocity dispersion M/R a a 2 , and galaxy 
structure (e.g., Strateva et al. 2001; Kauffmann et al. 2003; 
Blanton et al. 2003; Franx et al. 2008; Bell 2008). This corre- 
lation between the structural properties of galaxies with their 
stellar populations is important: it signals that the processes 
that determine the structures of galaxies at least correlate, and 
perhaps are the same as, the processes that shape whether 
or not a galaxy has cold gas and star formation. Further- 
more, these two populations evolve in their relative impor- 
tance: while the star-forming population has a stellar mass 
function that evolves slowly (Borch et al. 2006; Bell et al. 
2007; Peng etal. 2010; Brammeretal. 2011), the quiescent 
galaxy stellar mass function evolves rapidly from z ~ 2 to 
the present day (largely in normalization by factors of ~ 10, 
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but with modest or no evolution in shape or ‘characteris- 
tic’ mass AT; Bell et al. 2004; Borch et al. 2006; Faber et al. 
2007; Brown et al. 2007; Taylor et al. 2009; Peng et al. 2010; 
Dominguez Sanchez et al. 2011; Brammer et al. 201 1). 

A great deal of work, both theoretical and observational, 
has been carried out to try to better understand the drivers 
of the evolution of these two populations, particularly why 
some galaxies appear capable of shutting off their star for- 
mation while others are incapable of doing so. In this study, 
we will focus on processes that can shut down star for- 
mation in galaxies which reside in the center of their own 
dark matter halo. The clear effects of gas removal/starvation 
in dense environments (see, e.g., van der Wei et al. 2010, 
Pengetal. 2010, Weinmann et al. 2010, Pengetal. 2011, 
Weinmann et al. 2011 for recent discussions of this issue 
using survey datasets) are only a minor contributor to the 
evolution of ‘cosmic averaged’ galaxy population, owing to 
the small number of galaxies inhabiting dense environments 
(Pengetal. 2010). Accordingly, we do not discuss the ef- 
fect of environment in this paper (see, e.g., Peng et al. 2010 
for a careful discussion of the effects of environment as a 
function of cosmic time). An important point is that mod- 
els that include the growth of the dark matter framework, 
gas cooling, star formation and stellar feedback alone fail to 
predict a widespread population of non-star-forming galax- 
ies (Benson et al. 2003; Cattaneo et al. 2006; Somerville et al. 
2008; Dave et al. 201 1); all galaxies are expected to accrete 
(or cool) cold gas and form stars. 

A number of possible mechanisms have been proposed to 
keep a galaxy in the center of its own halo free of a signifi- 
cant cold gas content. Noting the strong tendency of quiescent 
galaxies to have prominent (or dominant) bulge components, 
it has long been thought that merging plays an important 
role in determining their structure (e.g., Toomre & Toomre 
1972; Barnes & Hemquist 1992; Hemquist 1993; Naab et al. 
2006; Hoffman etal. 2010). The approximate equality of 
the merger rate and the quiescent galaxy formation rate 
(e.g., Hopkins etal. 2010; Robainaetal. 2010), the detailed 
kinematic structure of early-type galaxies (Naab et al. 2006; 
Hoffman etal. 2010), and the empirical association be- 
tween relatively younger stellar populations and substructure 
(tidal tails, shells, asymmetries, etc.) in quiescent galax- 
ies (Schweizer & Seitzer 1992; Tal et al. 2009; Gyory & Bell 
2010) are qualitatively (and to a certain extent quantitatively) 
consistent with this picture. 

Largely in this merger context, the possibility that feedback 
from accretion onto a supermassive black hole may drive gas 
out of galaxies (Kauffmann & Haehnelt 2000; Springel et al. 
2005) or keep gas around galaxies from cooling (Croton et al. 
2006; Somerville et al. 2008) has been explored enthusias- 
tically. A wide array of observations are, at least at face 
value, qualitatively consistent with such a picture: e.g., galax- 
ies with big bulges have big black holes (Magorrian et al. 
1998; Giiltekin et al. 2009), low-redshift galaxies without 
prominent bulges cannot shut off their star formation on 
their own (Bell 2008), rapid large-scale winds are observed 
around post-starburst galaxies and quasars (Tremonti et al. 
2007; Prochaska & Hennawi 2009), and the energy measured 
in AGN-inflated outflow cavities in the hot gas atmosphere 
of groups and clusters of galaxies is approximately consis- 
tent with the energy required to offset cooling (Best etal. 
2006; Fabian etal. 2006). There are a number of other 
possible mechanisms for shutting off star formation, how- 
ever. A few examples are: the heating of the gaseous 


halo through virialization of the gas content (Naab etal. 
2007; Dekel & Bimboim 2008; Khochfar & Ostriker 2007; 
Johansson et al. 2009), changes in the mode of gas accre- 
tion onto galaxies as a function of dark matter halo mass 
(Kercs et al. 2005; Dekel & Bimboim 2006; Cattaneo et al. 
2006; Bimboim et al. 2007), and the heating of large halos 
with energy from cosmic rays (Guo & Oh 2008). 

1.1. The goal of this paper 

Given the range of possible mechanisms for shutting off 
star formation on galactic scales for galaxies in the center 
of their own halos, gathering empirical insight into the prop- 
erties of non star-forming galaxies, as a function of cosmic 
epoch, can be helpful. It has been argued that the key param- 
eter that correlates with the paucity of star formation is stellar 
surface density (Kauffmann et al. 2003; Franx et al. 2008), or 
possibly velocity dispersion (or roughly equivalently as M/R; 
Franx et al. 2008). Yet, for a sample of low-redshift galaxies 
from the SDSS, Bell (2008) instead argues that S6rsic index 
correlates much better with a lack of star formation for galax- 
ies in the center of their own halo (that could not have been 
stripped of their gas by external influences), as non-star form- 
ing galaxies have high Sersic indices but a range of surface 
densities that overlap with star-forming galaxies. Such an in- 
vestigation has not been carried out at z > 1 , owing to a lack 
of large-scale near-IR HST imaging until recently (see, e.g., 
Krieketal. 2009, Cameron etal. 2010, and Szomoruetal. 
201 1 for early progress towards this goal at z~ 2; Wuyts et al. 
201 lb explores this in more depth). 

A not inconsiderable challenge in achieving this goal is the 
definition of what constitutes star- forming or quiescent galax- 
ies. At a given stellar mass, the SFR of galaxies depends 
strongly on redshift, evolving by a factor of 5-10 by z~ 1 
and another factor of 4 or so out to z ^ 2 (Zheng et al. 2007; 
Noeskeetal. 2007; Dunne etal. 2009; Karim etal. 2011). 
Noting that the scatter in SFR at a given stellar mass is mod- 
est (a factor of < 2; Noeske et al. 2007), one could choose to 
define a quiescent galaxy as one that has a SFR more than la 
below the star forming galaxy locus at the redshift of inter- 
est, and a star-forming galaxy as any galaxy forming stars at 
a higher rate. An alternative approach is to separate galaxies 
by their optical-near- IR colors (e.g., U- VIV-J ; Wuyts et al. 
2007, Williams et al. 2009), where galaxies dominated by 
old stellar populations are distinguishable from star-forming 
galaxies with even substantial dust reddening (as used by e.g., 
Williams et al. 2009 and Brammer et al. 2011). In this paper, 
we adopt both techniques. We note that a galaxy at z ~ 2 
which is defined as quiescent according to these two criteria 
may have a SFR considerably in excess of almost all star- 
forming disk galaxies at the present day. While this means 
that our sample of quiescent galaxies does not have identi- 
cal properties across all redshifts, it does isolate a sample of 
galaxies with unusually low SFRs at that epoch given their 
stellar masses (confirmed by 24pm stacking) — one would 
like to understand why their SFRs are unusually low at that 
epoch. 

In this paper, we explore the evolution of the structures 
of galaxies as a function of redshift, and how they relate 
to the star formation activity in those galaxies. We use 
new near-infrared imaging from the Wide Field Camera 3 
(WFC3) on the Hubble Space Telescope taken as part of the 
Cosmic Assembly Near-IR Deep Extragalactic Legacy Sur- 
vey (CANDELS) multi-cycle treasury program (Grogin et al. 
2011; Koekemoer et al. 2011), focusing on the 0.6 < z < 
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2.2 galaxy population in the UKERT IR Deep Sky Survey 
(UKIDDS; Lawrence et al. 2007) Ultradeep Survey (UDS; 
§2). These data allow us to explore the structure of the 
M* > 3 >< 1O 1O M 0 galaxy population in the rest- frame optical 
to z< 2.2. We supplement this with data from the Sloan Digi- 
tal Sky Survey Data Release 2 (Abazajian et al. 2004) to con- 
nect with the properties of local galaxies. We perform two ba- 
sic analyses to explore the evolution of the galaxy population. 
First, we explore the evolution of the galaxy population drawn 
from an ‘equivalent’ constant comoving volume as a function 
of redshift, to get a sense of how the star formation and struc- 
tural properties of the galaxy population evolves with cosmic 
epoch (§3). Second, we explore the relationship between the 
structural parameters of galaxies and their star formation ac- 
tivity using the full sample at each epoch to maximize num- 
ber statistics, in an attempt to understand which structural pa- 
rameters best correlate with a lack of star formation activity 
(§4).In what follows, we use Vega magnitudes for rest-frame 
colors, assume that every star ever formed does so according 
to a universally-applicable Chabrier (2003) stellar IMF, and 
assume Hq = 70kms _1 Mpc' 1 , fl m>0 = 0.3 and 0 A ,o = 0.7. 

2. DATA 

2.1. UDS imaging data 

Our sample definition and selection is based on the 
public /f-band selected photometric catalog produced by 
Williams et al. (2009) 1 . We adopt these redshifts in this paper 
as the basis for conversion of apparent magnitudes into rest- 
frame colors, magnitudes and stellar masses, and the conver- 
sion of apparent to physical sizes. The Williams et al. (2009) 
catalog uses J and iC-band data for the UKIDDS UDS Data 
Release 1 (Lawrence et al. 2007; Warren et al. 2007) in con- 
junction with B, R, i and z-band imaging from the SXDS 
(Furusawa et al. 2008) and 3.6/nn and 4.5/im IRAC imaging 
data taken as part of the SWIRE survey (Lonsdale et al. 2003). 
Total fluxes were calculated using an elliptical Kron (1980) 
aperture, and observed-frame fluxes were calculated using 
a matched 1"75 circular aperture on PSF-matched images 
(with the obvious exception of the poorer resolution IRAC 
imaging data, whose fluxes were measured in 3" apertures 
and aperture corrected to the smaller aperture size). Finally, 
Williams et al. (2009) used the EAZY photometric redshift 
code (Brammer et al. 2008) to estimate photometric redshift 
from die photometric catalogs. 

Spectroscopic redshifts for the UDS are relatively few in 
number and preferentially focus on brighter sources; for these 
sources with ^ pec < 1.1 the redshift normalized median ab- 
solute deviation (the median absolute deviation, renormalized 
to give the same value as the RMS of a Gaussian distribu- 
tion) of |^>hot — -Zsp©c|/(1 + 3spec) is ~ 0.033 with 8% catas- 
trophic outliers (Williams et al. 2009). We have also com- 
pared the Williams et al. (2009) photometric redshifts against 
those of S. Wuyts et al. (in prep.), who used completely 
different (deeper) photometry and a similar photoz code to 
estimate photoz for galaxies in the CANDELS/UDS cover- 
age, finding a Az/(1 +zwiiii ams ) ~ 0.055, and 10% catas- 
trophic outliers (defined as having | Az|/(1 +zw imams) > 0.2. 
We have confirmed that use of the photometric redshifts, k- 

1 We choose to use this public catalog instead of other proprietary catalogs 

to better facilitate comparison with previous works and to allow easier repro- 

duction of the results. We confirm that the results and conclusions do not 
significantly change if repeated with the redshifts, colors and stellar masses 
from the currently proprietary catalog of S. Wuyts et al. (in preparation). 


corrections and stellar masses from Wuyts et al. instead of the 
public Williams et al. (2009) photometric redshifts plus the 
rest-frame colors and masses reported here yields no signifi- 
cant changes to our results. 

To explore the structure and morphology of the 0.6 < z< 
2.2 galaxy population, we use near-infrared F160W imag- 
ing from the HST, using the WFC3. CANDELS is a HST 
multi-cycle treasury program (Pis: S. Faber & H. Ferguson, 
PID: GO- 12060) to image five fields on the sky using the 
WFC3 and Advanced Camera for Surveys (ACS). The CAN- 
DELS imaging of the UDS field imaging includes 2/3 orbit 
in F125W and 4/3 orbits in F160W split into two epochs (see 
Koekemoer et al. 2011 and Groginetal. 2011 for more de- 
tails). A total of 44 WFC3 tiles were imaged in the UDS; 
when cross-matched with the Williams et al. (2009) catalog a 
total of 0.056 square degrees were covered. 

We use also lAjjm flux as a diagnostic to aid in the sepa- 
ration between galaxies with active star formation, and those 
with little or no star formation. We use 24/mi public data 
from the SpUDS survey 2 of the UDS. PSF fitting photometry 
with a 13" radius was performed, aperture corrected to total 
flux; the limiting flux of the catalog is ~ 50/Jy at the ~ 80% 
completeness (or 4cr level), and uncertainties in 24^m flux are 
of order ~ 20% (largely reflecting uncertainties in converting 
aperture to total flux, and source confusion). 

2.2. Rest-frame quantities and stellar masses 

The rest-frame magnitudes and stellar masses used in this 
paper were calculated using a set of template spectral en- 
ergy distributions from the PEGASE stellar population mod- 
els (see Fioc & Rocca-Volmerange 1997 for description of 
an earlier version of this stellar population model). Such 
stellar population models give very similar results to those 
of Bruzual & Chariot (2003), and stellar masses ~ O.lSdex 
higher than the models by Maraston (2005) 3 . These stellar 
population templates have solar metallicity, and we allow for 
dust attenuation following Calzetti (2001), with values of gas 
E{B- V) between -0.05 and 1.5 (to provide a little flexibil- 
ity to fit negative attenuation to account for small photometry 
problems, etc.; in practice small negative attenuation values 
are rare in the fits presented here). The templates included 
a broad range of exponentially-decreasing, constant or expo- 
nentially rising star formation histories, beginning at zf ~ 4 
(see Maraston et al. 2010 for a discussion of the importance 
of exponentially rising star formation histories for fitting z > 2 
SEDs). 

The templates treat the evolving galaxy population self- 
consistently, in the sense that all the galaxies that are in the 
z = 2 template set also appear in, e.g., the z= 1 or z= 0.6 
template sets further along their evolutionary path. Thus, 
the galaxy population at lower redshifts is required to have 
substantial older stellar populations, driving up somewhat the 
typical age of the stars and the typical stellar mass-to-light ra- 
tio at a given color. This is to be contrasted with other codes 
(e.g., Wuyts et al. 2008 or Pannella et al. 2009) which have 
redshift-independent template sets, but exclude those tem- 

2 http ://irsa. ipac.caltech.edu/data/SPITZER/SpUDS 

3 Although note that Krieketal. (2010) demonstrate that the 
Bruzual & Chariot (2003) models appear to fit the optical-near-ER SEDs of 
galaxies with large intermediate-age stellar populations substantially better 
than Maraston (2005) models do. The overall mass scale offsets between 
Bruzual & Chariot (2003) and Maraston (2005) models have been discussed 
by a number of papers, and recently by e.g., Brammer etal. (2011) and 
Domfnguez S&ichez et al. (2011). 
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plates that start star formation before the Big Bang — the typ- 
ical stellar populations of galaxies with this type of template 
fit end up being significantly younger than in the method dis- 
cussed here. It is unclear to us at this stage which method is 
more realistic; we will explore this issue with star formation 
histories drawn from a semi-analytic model of galaxy forma- 
tion in a future paper (see, e.g., Lee et al. 2010 for a version 
of this exercise at higher redshift). 

These templates were compared with photometric data 
points of each galaxy, given the photometric redshift value 
from Williams et al. (2009), and the template with the small- 
est x 2 was used to calculate rest-frame colors and a stellar 
M/L ratio. Rest-frame magnitudes were calculated by using 
the SED shape of the template to work out the predicted rest- 
frame magnitude of the object given the two nearest observed- 
frame bands, and then performing a weighted average of those 
two estimates of rest-frame magnitude. Stellar masses were 
estimated using the stellar M/Ls of the best-fitting template, 
referenced to the three longest observed wavelengths for over- 
all normalization. We note that the inclusion of dust, while 
it improves the quality of SED fit, leads to modest changes 
in the rest-frame colors and stellar masses of the sample ex- 
plored in this paper, given the overall degeneracy between 
the effects of dust extinction and stellar population age in the 
optical-near-IR spectral region; Bell & de Jong 2001). 

This method, when applied to galaxies with independently- 
estimated stellar masses (from independent photometry of 
similar but not identical datasets) in this field (S. Wuyts 
et al., in prep.) either the GOODS-S field (Wuyts etal. 
2008) or COSMOS field (Pannella et al. 2009), yields sim- 
ilar masses for intensely star-forming galaxies, and masses 
~ 0.2 dex larger for more dusty star-forming galaxies or non- 
star-forming galaxies (as the templates used here are more 
dominated by older stellar populations than the those used by 
Wuyts et al. 2008 or Pannella et al. 2009), with a scatter of 
0.2 dex. Rest- frame colors are reproduced to within 0.1 mag. 
When this method is compared with the masses of Bell et al. 
(2003). calculated on identical photometry and using a more 
restricted set of stellar population models without dust, there 
is no offset and 0.07 dex scatter in stellar masses, and a scatter 
of less than 0.05 mag in rest-frame U-V color. 

We adopt the stellar mass estimates described above for 
the purposes of this paper to ensure consistency of the stel- 
lar mass estimation method and stellar mass scale as a func- 
tion of redshift (the same code was used to estimate stellar 
masses at all redshifts, and the choice of templates evolves 
consistently from redshift to redshift). We have confirmed 
that the systematic discrepancies in stellar mass between the 
masses adopted here and those by, e.g., Wuyts et al. (2008) 
would operate to strengthen our conclusions (or in the case 
of the second part of the paper, leave them unaffected); their 
stellar masses for z > 0.6 galaxies are systematically lower 
by up to a factor of 2, with a factor of 2 scatter, and the 
evolution of the population would appear more rapid than 
it appears in this paper. In what follows we adopt a mass 
limit of 3 x 10 10 M 0 ; the sample is ‘complete’ above this 
limit 4 . The masses and rest-frame magnitudes of the sam- 
ple presented in this paper are available for download at: 

4 We 5 ay ‘complete’ because we filter the sample through steps that have 

uncertair:ies, e.g., photometric redshift or stellar mass estimation. We have 

confirmed that if we repeat the analysis with the Wuyts et al. sample, masses 

and rest-frame colors, drawn from deeper imaging data, our results and con- 

clusions do not significantly change. 


http://www.astro.lsa.umich.edu/~ericbell/data.php 

2.3. S&rsic profile fits of 0.6 < z < 2.2 galaxies 

To describe the structure of the galaxy population at 0.6 < 
z< 2.2, we use parametric S6rsic (1968) fits to the galaxy im- 
ages (A. van der Wei et al., in preparation). A surface bright- 
ness profile of the form E(r) = E e exp[-K,(jr l ^ n - 1 )] is fit using 
the GALFIT package (Peng et al. 2002), and the GALAPA- 
GOS wrapper 5 * * * , allowing the magnitude, axis ratio b/a , po- 
sition angle, half-light radius r e , Sersic index n, and central 
position to be free parameters. The Sersic index n describes 
the shape of the light profile, where n - 1 corresponds to an 
exponential light profile and n = 4 corresponds to a r 1 / 4 law 
profile characteristic of massive, spheroid-dominated early- 
type galaxies. The Sersic parameter n is a reasonably good 
proxy for the ratio of bulge luminosity to total luminosity, 
as illustrated in Fig. 14 of Simardetal. (2011) — systems 
with high n invariably host a prominent bulge, whereas sys- 
tems with low n host a weak or no bulge component. At the 
depths typical of this imaging, uncertainties in the fit parame- 
ters are 5\og l0 n ~ 0.15 dex, 5r e ~ 18%, and 8b/ a ~ 0.07, as 
constrained from independent GALFIT analyses on F105W 
imaging of a subsample of z ~ 1.6 galaxies in the UDS 
(Papovich et al. 201 1). 

The Sersic fits adopted in this paper are carried out on 
the F160W imaging data of CANDELS. This corresponds 
to rest- frame wavelength ranges of A rest ~ 0.55/0.65/0. 9/xm 
for z~ 2/ 1.4/0. 8 systems. A possible concern is that this 
change in rest-frame wavelength may affect the demograph- 
ics of the population. While file data at shorter wavelengths 
in the UDS field have not been analyzed with GALFIT, it is 
possible to test if this may be an issue using GALFITs on a 
smaller set of CANDELS F160W data in the GOODS South 
area, in comparison with published GALFITs on the F850LP 
ACS GEMS data for the extended Chandra Deep Field South 
(Haussler et al. 2007). We choose galaxies with 0.4 < z < 0.8 
for this test, where the F850LP data span the same range in 
rest-ffame wavelength as the F160W data for 1.3 < z< 2.2. 
For systems with low n < 1, we find a slight tendency for 
the F160W Sersic index to exceed the F850LP data (by <0.1 
dex), and for n > 2 there is no systematic difference between 
the two sets of Sersic fits. The scatter around these modest 
offsets is ~0.2 dex, equivalent to the combined uncertainties 
of the fits. The fraction of systems with n > 2.5 in F850LP 
(revStframe ~ 0.6/zm) is in fact 20% larger than the fraction de- 
rived using the F160W imaging (restframe ~ 1/im), opposite 
in sign and much smaller than the factor of several observed 
change in the number of high-n galaxies from z ~ 2 to the 
present day. We conclude that our use of F160W data alone 
across the 0.6 < z < 2.2 redshift range is not an important 
source of systematic error in this analysis. 

2.4. SDSS parameters for the low-redshift comparison 
sample 

In order to connect with the present-day galaxy population, 
we use a sample of low-redshift galaxies explored in Bell 
(2008) from the SDSS Data Release Two (Abazajian et al. 
2004), and presented in the NYU Value-Added Galaxy Cat- 
alog (Blanton et al. 2005). We use foreground extinction- 
corrected (Schlegel et al. 1998), /r-corrected (Bell et al. 2003) 

5 See http://astro-staff.uibk.ac.at/~m.barden/galapagos/ 
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r-band absolute Petrosian magnitude for the galaxy ab- 
solute magnitude (random and systematic uncertainties < 
0.15 mag), and model colors for higher S/N estimates of 
galaxy color (uncertainties < 0.05 mag). Following Bell et al. 
(2003), we have merged this catalog with the Two Micron 
All-Sky Survey (2MASS; Skrutskie et al. 2006), to facilitate 
SED fitting and to allow splitting of galaxies by U- VJV-J 
into quiescent and star-forming populations. Stellar masses 
and rest- frame colors were estimated from ugrizJK photome- 
try using the above stellar population model templates. 

Star formation and AGN classifications and estimates of 
total star formation rate were taken from Brinchmann et al. 
(2004), using emission line measurements described in 
Tremonti et al. (2004). Galaxies are classified as star forming, 
AGN, composites, or are left unclassified (typically because 
the galaxies lack line emission in their SDSS spectra). 

The only source of Sersic fits for all galaxies in our sam- 
ple is Blanton et al. (2003) 6 * * * . Blanton et al. (2003) fit the light 
profile of galaxies in the SDSS, measured in circular aper- 
tures, with a seeing-convolved Sfrsic (1968) profile for all of 
the galaxies in the VAGC. The S6rsic fits by Blanton et al. 
(2003) give values for the Sersic index, in particular, that 
are offset from n values determined using 2D galaxy im- 
age fits. We have compared the S6rsic indices (and other 
fit parameters) from Blanton et al. (2003) to fits carried out 
by van der Wei (2008) on a small subset of galaxies in the 
NYU VAGC. We find that the S6rsic index estimates are 
related: log 10 /72D ~ -0.39 + I.751og 10 /7NYUVAGC» with 0.2 
dex scatter. Half-light radii show the following correlation: 
log 10 ^D ~ log | 0 /Kyu VAGC “0.05 +0.025i72D,new> with 0.1 dex 
scatter, where /? 2 D,new is the estimate of equivalent 2D Ser- 
sic index derived from the NYU VAGC S6rsic index using 
the above relation. A similar analysis was carried out with 
completely independent GALFITs by Guo et al. (2009), and 
importantly the above trends are identical in the case of n, 
and < 0.05 dex different in the case of r, to the median offsets 
as a function of n in their Fig. Al. Recall that the 0.2 dex 
scatter between the ‘rescaled’ NYU n values and those of 
van der Wei (2008) or Guo et al. (2009) (or the 0.1 dex scatter 
in radii) is comparable to the typical joint uncertainties in any 
comparison of even unbiased values of n (or r). We conclude 
that these rescaled NYU n and r values are unbiased, have un- 
certainties comparable to those determined directly from 2D 
fits, and are appropriate for connecting the evolution of n and 
r with the results of 2D fitting for the z > 0.6 galaxies. 

3. THE EVOLVING RELATIONSHIPS BETWEEN STELLAR MASS, 
COLOR AND MORPHOLOGY 

3. 1 . Separating galaxies into quiescent and star forming 
populations 

We choose to separate galaxies by their star formation ac- 
tivity using two sets of independent diagnostics: mid infrared 
(or for the SDSS emission line) information, and position on 
optical-near-infrared color-color diagrams. We use both cuts 
in this paper; we describe the 24/im-derived cuts first (show- 
ing how they relate to rest-frame color cuts) and then show 
the rest-frame color cuts (showing how they relate to 24pm 

6 Simard et al. (201 1) also fit all galaxies in the SDSS with single Ser- 

sic profile fits. We find that for our particular sample half of the fits (those 

with larger quoted formal uncertainties from Simard et al. 2011, we believe 

indicative of failed fits) have very large discrepancies compared to fits from 
Blanton et al. (2003) or van der Wei (2008). Because these poor fits comprise 

some ~ 1/2 of our sample, we cannot adopt the estimates of Simard et al. 
(2011) for this work. 


cuts). 

At z>0.6, galaxies are classed in large part according to 
their 24/im emission properties as an admittedly imperfect 
proxy for obscured SFR (we use emission-line diagnostics 
and SFR estimates from Brinchmann et al. 2004 for galaxies 
from the SDSS). Substantial 24/zm flux can also result from 
AGN activity. We do not attempt to discriminate between 
AGN activity and SF activity for our purposes here, simply 
noting that systems with their 24pm flux dominated by AGN 
at this 24pm luminosity and redshift range are not the dom- 
inant population (e.g., Donley etal. 2008; Kartaltepe et al. 

2010) , and that we are primarily attempting to weed out galax- 
ies whose rest- frame optical colors are a poor reflection of the 
stellar populations in that galaxy, a goal for which our simple 
approach is sufficient. 

In this spirit, we wish to avoid an explicit, and uncertain, 
conversion of 24pm flux into SFR (e.g., Papovich et al. 2007; 
Elbaz et al. 2011; Wuyts et al. 201 la). At each redshift of in- 
terest, we fit the relationship between 24pm flux and stellar 
mass for galaxies detected at 24pm (shown as filled symbols), 
as shown in Fig. 1. The approximate trend at all redshifts is 
i°gio W/*Jy *V/ 2 + 0.51og 10 (M*/3 x 10 10 Afo), with a scat- 
ter of less than 0.3 dex (we use the actual fits and scatters, 
which vary slightly with redshift, to perform the split into star 
forming and quiescent). The slope and scatter of this relation- 
ship is well-documented and studied (e.g., Salim et al. 2005; 
Zheng et al. 2007; Noeske et al. 2007; Karim et al. 201 1); it 
is a remarkable coincidence that the zero point in terms of 
24pm flux varies so little with redshift, owing to the inter- 
play between the dramatic reduction of SFR at a given stel- 
lar mass with decreasing redshift, the luminosity distance and 
the redshift-dependent 24pm ^-correction. Quiescent galax- 
ies then must have a 24pm flux (UDS) or SFR (SDSS) lower 
than —1(7 from the star forming galaxy locus. 

In Fig. 1, we have color-coded symbols by their position on 
the rest-frame U- VIV-J diagram (Fig. 2), using the slightly 
redshift-dependent cuts described in Williams et al. (2009). 
Orange symbols show galaxies with rest-frame optical-near- 
IR colors characteristic of quiescent galaxies, and purple sym- 
bols show galaxies with colors characteristic of star-forming 
galaxies with a range of reddening values. In Fig. 2, we show 
the optical-near-IR colors of galaxies in the six redshift inter- 
vals of interest in this paper. In contrast to Fig. I, we have 
color-coded the symbols in Fig. 2 by 24pm flux (UDS) or 
SFR (SDSS). Galaxies with 24pm fluxes/SFRs lower than 
-1 <j from the star forming galaxy locus have been color- 
coded red, and galaxies with fluxes/SFRs higher than -1 o 
from the SF galaxy locus have been color-coded blue. In the 
z < 0.05 slice, emission-line diagnostics are available, and 
any object with AGN- like lines or composite SF/AGN lines 
(Brinchmann et al. 2004) have been color-coded green. 

Inspection of Figs. 1 and 2 shows the large degree of 
overlap, and the complementarity of having both explicit 
24/im/SFR information and U- VIV-J colors (see also 
Williams et al. 2009, Wuyts etal. 2009 and Brammeretal. 

2011) . Galaxies with quiescent U- V1V- J tend, for the 
most part, to be undetected at 24pm (Fig. 1). There are ex- 
ceptions to this: a few galaxies with quiescent U- VIV-J 
at z < 1.5 are detected at 24/im, and by z ~ 2 it is clear 
that for our particular dataset the contamination of the qui- 
escent region of U- VIV-J color space by 24/im-detected 
objects is significant. Stacking at 24pm of the remaining 
individually-undetected galaxies with quiescent U- VIV-J 
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FIG. I. — Flux at 24/xm as a function of stellar mass in three different redshift bins; detections are shown as filled symbols and 4a upper limits shown by 
arrows. The thick solid line and dotted lines are the robust linear fits to the data and the ±la range. Colors denote classification according to U- VIV-J (Fig. 
2): purple symbols denote galaxies that lie in the star forming part of the color-color diagram, and orange galaxies line in the quiescent region of the plot. The 
large solid purple points with error bars show the stacked 24/xm flux and uncertainty for galaxies individually undetected at 24/xm but in the star-forming part of 
the U- V/V-J diagram. 



0.5 1.0 1.5 2.0 2.5 0.5 1.0 1.5 2.0 2.5 

V-J (rest-frame) V-J (rest-frame) 


Fig. 2. — Rest-frame V- V as a function of V-J color for the six redshift 
bins used in this paper. Galaxies are color-coded according to their 24/im 
properties (z> 0.6) or emission line properties (z< 0.05); (quiescent) galax- 
ies not detected or fainter than -la from the star forming galaxies locus are 
color-coded red, and the rest of the galaxies (all star forming) are color coded 
blue. Green galaxies at z < 0.05 are galaxies in the SDSS with AGN-like 
emission lines. Superimposed are the rest-frame color cuts used in this paper, 
following Williams et al. (2009). 


yields marginal (2-3a significance) detections at the 5- 10/xJy 
level at all redshifts (indicating SF/AGN activity a factor of 
>10 lower than typical star-forming galaxies at that redshift; 


see also Papovich et al. 2006 for an early discussion of star 
formation in red-selected galaxies). 

Conversely, galaxies detected clearly at 24/xm are almost 
always in the star forming region of U- V/V-J (Fig. 2), but 
again with some exceptions (e.g., at z~ 1.4 there is a clear 
group of galaxies with star forming colors but are individually 
undetected at 24/xm). Stacks of those few individually 24/xm 
undetected galaxies with star-forming colors yields significant 
detections at the 15-25/zJy level (blue filled points with error 
bars on Fig. 1), a factor of a few below the SF galaxies locus, 
largely consistent with an interpretation of these systems as 
the low SFR tail of the SF galaxy population. This high degree 
of correspondence between the two methods has been shown 
before by e.g., Williams et al. (2009) and Wuyts et al. (2009). 
Fig. 2 also shows that the use of rest-frame color informa- 
tion for z < 0.05 is particularly valuable; galaxies classified 
as AGN can have either quiescent or star forming colors. 

We choose to separate galaxies into quiescent and star form- 
ing using both criteria, to capitalize on their different strengths 
and shortcomings. Discrimination by U- V/V-J is sensi- 
tive to lower amounts of star formation than 24/xm separa- 
tion, especially at z~ 1.4 and z~ 2. On the other hand, 
separation by 24/xm is considerably less sensitive to photoz 
error than U -V/V-J, as one simply needs to know which 
redshift bin the galaxy is in, and even some cross-talk be- 
tween bins can be tolerated. We define quiescent galaxies 
as having both ‘quiescent’ colors in U- V/V-J and 24/xm 
fluxes/SFRs lower than -la from the star forming galaxy lo- 
cus at the redshift of interest. For SDSS galaxies with AGN- 
like emission lines, we split only on the basis of U- V/V-J. 
Star forming galaxies are defined as those which satisfy ei- 
ther (or both) of the U- V/V-J star forming galaxy color 
cuts or having 24/xm fluxes/SFRs brighter than - la from the 
star forming galaxy locus. One can see that incorporating 
24/xm data and insight from U- V/V-J into this analysis is 
crucial. As can be seen directly in Figs. 1 and 2, and fur- 
ther appreciated by the intermixing of blue and red symbols 
at red rest-frame U- V in Figs. 3-6, failure to flag galaxies 
by multiwavelength-derived SF activity leads to considerable 
confusion between star forming galaxies with substantial dust 
columns and non star-forming galaxies, especially at higher 
redshifts (Taylor et al. 2009; Brammer et al. 201 1). 
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Fig. 3.— The evolution of U- V rest-frame color (in Vega magnitudes) as a function of stellar mass, in six different redshift bins. Galaxies that 
appear in this figure are subsampled to a volume of ltfMpc 3 ; i.e., this shows the evolution of the massive galaxy content of a ‘representative’ 
volume >vith cosmic time. Open symbols show galaxies with n < 2.5 and filled symbols show galaxies with n > 2.5. For all symbols the axis 
ratio of the galaxy is reflected by the axis ratio of the symbol, and the size of the symbol scales with (1 + log 10 r e / kpc). The black line is shown 
in all panels at the approximate locus of 2 red galaxies. Galaxies at z> 0.6 are color-coded by star formation activity — galaxies classified 
as quiescent by both 24um/SFR and U- VIV- J diagnostics are color-coded red, and blue symbols show galaxies classified as star forming 
using either u- VfV-f or 24jxm/SFR. AGN with z< 0.05 are classified using only U- VIV-J. In all panels, the small red arrow shows the 
median mass of the whole (not subsampled to lO^pc 3 ) sample at that redshift, and the small red numbers on the left/right show the fractions 
of quiescent galaxies in the halves of the whole sample with masses below/above that median mass respectively. 
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Fig. 4 — The evolution of U-V rest-frame color (in Vega magnitudes) as a function of 0.5M/r e ~ <r /2 , in six different redshift bins. Symbols are as in Fig. 
3 for z > 0.6. For the SDSS, we code galaxies only by their star formation activity for clarity. The first five panels (z > 0.6) are for the full UDS survey and 
the SDSS panels use the whole SDSS DR2 subsample used in Bell (2008). In all panels, the small red numbers on the left/right show the fractions of quiescent 
galaxies in the halves of the sample below and above the median a* (red arrow), respectively. 


3.2. Evolution of the galaxy population in a given comoving 

volume 

Fig: 3 shows the rest-frame U-V colors of galaxies, as a 
function of their stellar mass, in six different redshift bins. We 
choose to show the properties of the galaxy population as a 
function of U- V rest-frame color to connect with other stud- 
ies (e.g.. Bell et al. 2004; Borch et al. 2006; Ruhland et al. 
2009; Whitaker et al. 2010; Brammeretal. 2011), and as a 
joint (rather sensitive) constraint on SFH and dust content. 
Galaxies are color-coded by star formation activity (§3.1): 


red symbols show galaxies classified as quiescent using both 
24pm information and U- V/V- J colors, and blue symbols 
show the star forming galaxy population. In all panels of this 
figure, the galaxy population has been Monte Carlo subsam- 
pled down to an equivalent volume of 10 5 Mpc 3 by adjusting 
the number of galaxies to track the number density of galax- 
ies with M* > 3 x 10 10 M g determined from larger surveys 
(the line in Fig. 1 1). Put differently, variations in the num- 
ber of galaxies from panel to panel illustrate true evolution in 
the galaxy population (as the volume is fixed; see Appendix 
A for further discussion). Filled symbols show galaxies with 
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samples and symbols are as in Fig. 4. In all panels, the small red numbers on the leftAight show the fractions of quiescent galaxies in the halves of the sample 
below ar.d above the median stellar surface density (red arrow), respectively. 


n > 2.5. Open symbols show galaxies with n < 2.5. In all pan- 
els, the linear size of the symbol scales with (1 +log 10 r e /kpc), 
and the axis ratio of the symbol is the same as that of the 
galaxy of interest. The black line is shown in all panels, 
for reference, at the approximate position of z ~ 2 non star- 
forming galaxies. 

The evolution of the galaxy population in the epoch z~ 2 
to the present day is obvious. As has been argued by a num- 
ber of other authors (e.g., Amouts et al. 2007; Fontana et al. 
2009; Taylor etal. 2009; Ilbertetal. 2010; Cassataetal. 
2011; Dommguez Sdnchez et al. 201 1 ; Brammer et al. 201 1), 
there is dramatic evolution in the overall number of galaxies 


with M* > 3 x 10 10 M o (as quantified in Fig. 11). Further- 
more, Fig. 3 shows that the evolution of the number density of 
quiescent galaxies is particularly striking (again, as has been 
argued by the above cited works). Fig. 3 makes it clear how- 
ever that the evolution of the star formation activity of the 
intermediate-mass and massive galaxy population is accom- 
panied by a large-scale change in the structure of galaxies 
(see also Wuyts et al. 201 lb, who see similar behavior) from 
a z> 1.5 population dominated by low n (little or no bulge), 
mostly star-forming systems to the present population, domi- 
nated by galaxies with high n (with a prominent bulge), many 
of them quiescent (but not all of them). Quantitatively, there 
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Fig. 6. — The evolution of U-V rest-frame color (in Vega magnitudes) as a function of S6rsic index, in six different redshift bins. The samples and symbols 
are as in Fig. 4, except that now filled/open symbols denote those galaxies with above/below median surface densities at that redshift (in the stellar mass range 
considered in this paper). In all panels, the small red numbers on the left/right show the fractions of quiescent galaxies in the halves of the sample below and 
above the median Sersic index (red arrow), respectively. 


are a factor of ~2.5 (3.5) more quiescent ( n > 2.5) galaxies 
with M* > 3 x 10 10 Mq today than there were galaxies with 
those masses at z~ 2, respectively. Thinking about it differ- 
ently, the current population of quiescent ( n > 2.5) galaxies 
with M* > 3 x 10 1 o A/q is approximately as numerous as the 
entire Af* > 3 x 10 ,o M© population at z~ 1.1 (0.7) respec- 
tively. This change in global demographics from 2 to the 
present day makes it clear that in addition to processes that 
shut off star formation on galactic scales, there must be (the 
same or different) processes that lead to an associated change 
in the surface brightness profiles of galaxies over the same 


time period (and given the good correspondence between a 
lack of star formation and structure, arguably with similar 
time scales). 

Fig. 3 shows also that the scatter in the quiescent galaxy 
CMR decreases towards lower redshift TTie evolution of 
CMR scatter from z~ 2 to the present day is well-documented 
in the literature (Ruhland et al. 2009; Whitaker et al. 2010). 
The scatter in the U-V colors of non star-forming galax- 
ies with U- V > 0.6 at z> 1.3 is 0.17 mag (our measure- 
ment), very consistent with the carefully-measured results of 
Whitaker et al. (2010), who find a scatter of 0.13-0.2 mag for 
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1 .3 < 2 < 2. Ruhland et al. (2009) find that the scatter in U- V 
color at z < 1 is ~ 0. 1 mag (measured much more carefully 
than the CMRs presented in Fig. 3; our measurements also 
give a scatter of ~ 0.1 mag), essentially independent of red- 
shift. Modeling presented in both Ruhland et al. (2009, for 
z< 1) and Whitaker et al. (2010, at 1 < z < 2) shows that 
the evolution of CMR scatter is naturally interpreted as be- 
ing caused by a constant inflow of new galaxies on to the red 
sequence at the observed number density growth rate. 

4. EMPIRICAL CORRELATIONS BETWEEN A LACK OF STAR 
FORMATION AND GALAXY STRUCTURE 

4.1. Broad trends 

One of the most notable trends seen by Amoutsetal. 
(2007), Taylor etal. (2009), Ilbertetal. (2010), 
Dominguez Sanchez et al. (20 1 1), Brammer et al. (20 1 1) 
and elaborated upon in Fig. 3 is the dramatic growth of the 
quiescent galaxy population from z ~ 2 to the present day. 
Fig. 3 demonstrates also that with the growth of the quiescent 
galaxy population comes a concurrent growth of the popula- 
tion of concentrated n > 2.5 galaxies (see also Wuyts et al. 
2011b). Franxetal. (2008) argues that these galaxies also 
have high surface density and ‘velocity dispersion’, and that 
‘velocity dispersion’ is the parameter that best correlates with 
a lack of SF activity at z < 2. In this section, we explore 
how the different structural parameters correlate with star 
formation activity in an attempt to gain possible insight into 
the processes that drive galaxies into quiescence. 

In Figs. 4 and 5, we show the U- V rest-frame colors of 
galaxies as a function of 0.5M/r e oc cr /2 (following Franx et al. 
2008, where cr f is used to differentiate this quantity from a 
true velocity dispersion; it should be noted that this quan- 
tity also scales with rotation velocity, with a different pro- 
portionality constant, for systems dominated by rotation), and 
0.5JW/tj^ = E (surface density within the half-light radius, 
following Kauffimann et al. 2003, and as explored for low red- 
shift by Bell 2008). Symbols are coded by SSrsic index, axis 
ratio, size and star formation activity as in Fig. 3. We draw 
from the full UDS and SDSS DR2 datasets, as opposed to 
sub-sampling down to an equivalent volume of 10 5 Mpc 3 , in 
order to delineate the trends with better fidelity than the sub- 
sample shown in Fig. 3. 

In Figs. 4 and 5, one can see trends reported by 
Kauffmannetal. (2003) and Franxetal. (2008): galaxies 
with either high velocity dispersion, or high surface density, 
tend not to form stars. Yet, in these figures one can also see the 
behavior reported by Bell (2008): there are a number of non- 
star forming galaxies also at low mass, velocity dispersion or 
surface density, and to a great extent their main defining char- 
acteristic is a relatively high S6rsic index 7 * * . 

Fig. 6 shows the trend in rest-frame color with Sersic index. 
This trend was not explored by Franx et al. (2008), as they did 
not analyze large-scale HST imaging, and therefore lacked re- 
liable measurements of surface brightness profile shape; this 
can, however, be compared with (and is consistent with) Fig. 
1 of Wuyts et al. (2011b). Symbols are similar to previous 
figures, except that filled/open symbols now denote galaxies 

7 The actual fraction of galaxies forming stars at low S or cf may be rather 

higher, as the sample is limited by stellar mass. Fig. 3 shows that quiescent 

galaxies are smaller at a given stellar mass than star-forming galaxies, there- 

fore it is possible that if Figs, 4 and 5 were ^ or E limited samples they 
would show a more prominent population of (lower stellar mass) quiescent 
galaxies with relatively low er' or 


with above/below median stellar surface density at that red- 
shift for the galaxies with A/* > IO lo * 5 M 0 in our sample (re- 
call that the previous filled/open distinction was by Sersic in- 
dex, which would be redundant). Focusing on the z > 0.6 
points, one can see that galaxies with high Sersic index are 
much more likely to be non star-forming than their low Ser- 
sic index counterparts. Furthermore, one can see that there 
is a range of stellar surface densities at a given S6rsic index. 
A similar qualitative behavior is seen for the SDSS galaxies, 
notwithstanding quantitative differences in the definition of 
S6rsic index for the NYU VAGC catalog. Put differently, the 
degree of left-right split on the plots between red and blue 
points (lack of detectable star formation vs. star formation) 
appears stronger in Fig. 6 than in Figs. 3-5. 

4.2. Which parameter correlates best with a lack of star 
formation? 

Figs. 3-6 demonstrate that ‘typical’ quiescent galaxies have 
higher mass, ‘velocity dispersion’, surface density and Sersic 
index than ‘typical’ star-forming galaxies. In this section, we 
will explore further which parameter correlates the best with 
a lack of star formation activity. 

Fig. 7 shows the quiescent fraction (evaluated in miming 
bins of 101 galaxies at z> 0.6, or 501 galaxies at z< 0.05) as 
a function of the rank of a galaxy in stellar mass (dotted line), 
‘velocity dispersion’ (grey solid line), surface density (dashed 
line) and S6rsic index (black solid line) in three broad redshift 
bins. At all redshifts, stellar mass is a poor predictor of qui- 
escence. At z < 0.05, Sersic index is clearly a better predictor 
of quiescence than any other parameter; in particular, galaxies 
with low n overwhelmingly host detectable star formation. 

At 0.6 < z < 2.2, one can see a slightly different situation: 
‘velocity dispersion’, surface density, and Sersic index corre- 
late comparably well with SF activity. Yet, in both bins, one 
can see a weaker version of the behavior seen at low redshift: 
S6rsic index correlates somewhat better with a lack of star 
formation than either ‘velocity dispersion’ or surface density. 
The small numbers on the lower left and right hand comers 
of each panel of Figs. 3-6 help to illustrate this point. These 
numbers show the fraction of quiescent galaxies for two dif- 
ferent subsamples: the half with lowest stellar M*, <r', E, or 
n y and the half with highest A/*, cr', E, or n. The S6rsic in- 
dex is the metric that tends maximizes the contrast between 
the two halves of the sample (except at 0.9 < z< 1.3 and 
1.5 < z < 1.8, where Sdrsic index still correlates very well 
with quiescence). Note that repeating this analysis with al- 
ternate photometric redshifts, stellar masses and rest-frame 
properties (S. Wuyts et al., in prep.) yields a slightly stronger 
preference still for n as the parameter that best correlates with 
quiescence. 

Fig. 8 emphasizes this finding. In the same three broad red- 
shift bins, we show U-V color of galaxies as a function of 
‘velocity dispersion’ (found by Franx et al. 2008 and E. Che- 
ung et al., in prep., to correlate marginally better with a lack 
of star formation than surface density). Symbols are coded 
as they are in Figs. 3-5. The top panels show galaxies with 
n< 1.5 (galaxies with little or no bulge) and the bottom panels 
show galaxies with n > 3 (galaxies with a prominent bulge). 
The population with n < 1.5 is overwhelmingly star fonning. 
In strong contrast, the n > 3 population has a large quiescent 
fraction. Strikingly, one can see that star formation or quies- 
cence happens over a wide range in ‘velocity dispersion’ . 

There are two emergent themes that we wish to draw the 
reader’s attention to. First, Figs. 6 and 8 show that, with very 
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Rank Rank Rank 

Fig. 7. — The quiescent fraction in three broad redshift bins as a function of the rank of n, E, <j\ and M+ (rank denoting where a galaxy is in the sorted list of 
the /j, o' or M* in the sample at that redshift of interest). These are running average fractions determined for the ±50 (250 for the SDSS) nearest neighbors in 
rank of the quantity in question. Mass is always a poor predictor of quiescence; Sersic index is clearly superior to E or <r' at z < 0.05, and performs as a slightly 
better predictor of quiescence for all z < 2.2. 


few exceptions, galaxies with low Sersic index all appear to 
form stars at all z < 2.2. The threshold appears to be some- 
where around n ~ 1.5-2: at n < 2, the fraction of quiescent 
galaxies is < 10% (and in many redshift bins it is less than 
a few percent). When investigated in more detail using Che 
same SDSS sample (Bell 2008), it was found that i) real low 
n quiescent galaxies are all satellite galaxies in galaxy clus- 
ters, i.e., they are stripped disk galaxies, and ii) the few qui- 
escent ‘low ri systems in the centers of their own halos that 
remained were in fact the result of measurement error in n, 
as visual inspection showed a distinct bulge component. We 
show examples of some n <2 quiescent galaxies at z> 1.5 
in Fig. 9 (at 0.9 < z < 1.5 there are only three n <2 qui- 
escent systems, and those look similar to the z> 1.5 exam- 
ples; the 0.6 < z < 0.9 points are bad S6rsic fits). While some 
systems are relatively extended and have low n, and the one 
inclined galaxy is clearly reminiscent of a disk, most appear 
spheroidal and compact. Given Figs. 6 and 9 in concert, it is 
clear the vast majority of quiescent galaxies have a prominent 
spheroid. This extends the results of Bell (2008) determined 
for nearby galaxies to 2.2, when the Universe was ~ 1/4 
of its present age: galaxies lacking a prominent bulge appear 
to have great difficulty shutting off their own star formation 
on galactic scales. 

Second, Figs. 6 and 8 make it clear that having a high Ser- 
sic index alone (or indeed, having high n, cr 7 , E and M*) is not 
enough to ensure a lack of star formation. At all redshifts, a 
small minority of high n sources form stars at an appreciable 
rate. This illustrates a key point of this paper; it appears that 
for all z < 2.2, a large bulge is necessary to stop star forma- 
tion, but is not sufficient to stop star formation. This extends 
the conclusion of Bell (2008) determined for local galaxies 
out to z ~ 2.2, when the non star-forming galaxy population 
was considerably less prominent. 

5. DISCUSSION 

There are two main observational results in this paper: the 
rapid growth of the quiescent galaxy population between z- 
2 and the present day, and the recognition that this growth 
appears to be intimately linked to the growth of the high S6rsic 


index galaxy population. 

5.1. Musings on the mechanisms tha t prevent significant cold 
gas in galaxies 

These results have some bearing on understanding which 
mechanisms lead to quenching of star formation in galaxies. 
Recall that the role of environmental quenching is relatively 
minor in our ‘cosmic averaged’ population evolution, and that 
we are focusing on which types of physical process lead to 
quenching of star formation in galaxies in the centers of then- 
halos (the mass quenching of Peng et al. 2010). To facilitate 
this, we will set up two straw person hypotheses: suppression 
of star formation by feedback (either star formation or AGN 
feedback; ‘feedback quenching’, e.g., Kauffmann & Haehnelt 
2000, Croton et al. 2006, Somerville et al. 2008), or suppres- 
sion of star formation because the halo reaches a certain crit- 
ical mass (‘halo quenching’; e.g., Dekel & Bimboim 2006, 
Cattaneo et al. 2006). 

Fig. 7 demonstrates that internal properties (Sersic index, 
and to a lesser extent ‘velocity dispersion’ and surface den- 
sity) correlate strongly with quiescence, whereas galaxy mass 
only weakly correlates with star formation activity. Galaxy 
mass is expected to correlate well with halo mass. More et al. 
(2009) measure the scatter in luminosity at a given halo mass 
to be 0.16dex, substantially less than the ~ ldex dynamic 
range probed in this work. The observed weakness of cor- 
relation between quiescent fraction with stellar mass, coupled 
with the expected modest scatter between stellar mass (via its 
proxy, luminosity) and halo mass, implies a weak correlation 
between halo mass and quiescence at all z< 2.2. This is con- 
sistent with the claim by More et al. (201 1) that at fixed stellar 
mass there is no difference between the average halo masses 
of quiescent and star-forming central galaxies. We conclude 
that quiescence is not determined by halo mass alone. 

The results here demonstrate that systems with high Ser- 
sic index, and to a certain extent systems with high ‘velocity 
dispersion’ and high surface density, are much more likely 
to be quiescent (see also Franx et al. 2008 and Wuyts et al. 
201 1 b). In the context of gas expulsion (or prevention of cool- 
ing) by feedback, it is suggestive that the best correlation with 
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Fig. 8. — The evolution of U- V rest- frame color (in Vega magnitudes) as a function of 0.5A//r e ~ cr' 2 , in three broad redshift bins. The top panels show 
galaxies with n < 1.5 (galaxies with little or no bulge), the bottom panels show only galaxies with n > 3 (galaxies with a prominent bulge component). Symbols 
are as in Fig. 4. A lack of star formation (red symbols) can happen over a wide range in a * , but correlates strongly with n. 


quiescence is SSrsic index (an imperfect proxy for the rela- 
tive prominence of a bulge component). Given the black hole 
mass-bulge mass correlation (Haring & Rix 2004; Peng et al. 
2006; Giiltekin et al. 2009), our results tentatively support the 
AGN feedback paradigm, at least at the qualitative level. In 
particular, the seeming inability of galaxies with low Sersic 
indices to shut off their star formation is very naturally in- 
terpreted in this ffamework (Bell 2008) — no supermassive 
black hole, no shut-off of star formation. We note that feed- 
back from star formation-driven winds may also be of rele- 
vance, but we caution that the results shown here argue that 
such winds must be efficient at wholesale ISM removal only 
during star formation events that create a prominent bulge (as 
disk galaxies keep forming stars with enthusiasm). 

5.2. Quiescent disks? 

There has been much recent discussion of the pres- 
ence and importance of stellar disks in quiescent galaxies. 
At z> 1.5, Stockton et al. (2004), McGrath etal. (2008), 
van Dckkum et al. (2008) and van der Wei et al. (201 1) have 
argued that most quiescent galaxies with masses in excess of 
10 n M o have prominent stellar disks. Bundy et al. (2010) has 
explored this issue at z < 1.2, finding that a large fraction of 
quiescent galaxies have disks in addition to significant bulge 
components; indeed, the existence of SOs and disky ellipticals 
in all environments is well-known (see, e.g., van den Bergh 
2009, and references therein). The highly-inclined fraction of 
these systems are clearly visible in Fig. 6, as the elongated 



FIG. 9. — F160W Postage stamps of quiescent galaxies with n < 2 and 
stellar surface densities between 10 9 and 10 10 A/q kpc“ 2 . The postage stamps, 
within each class, are ordered by stellar mass (ordered left to right). At all 
redshifts, all postage stamps are 40 physical kpc on a side, and are scaled 
to a constant ’stellar mass density’ (total intensity is scaled to total stellar 
mass, meaning that if the stellar M/L is constant over the face of the galaxy 
this postage stamp should reflect the stellar mass density) and are displayed 
using asinh scaling (linear at low intensity, and logarithmic at higher surface 
brightness; Lupton et al. 1999). 

symbols with low b/ a. Such systems appear to be somewhat 
more common at z > 1 .5, but are present (especially at lower 
Sersic indices n~ 2-3) at all redshifts. 

There are two comments that we wish to make about this 
issue. Firsdy, the vast majority of these systems have Sersic 
indices n > 2. In McGrath et al. (2008) and van der Wei et al. 
(2011), bulge/disk decompositions were carried out, and 
showed that these systems with relatively high Sersic in- 
dex are also well-explained as composite bulge/disk systems 
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with relatively large bulges. Furthermore, kinematic stud- 
ies of local quiescent galaxies has demonstrated that the 
vast majority of quiescent galaxies have significant rotation 
(Emsellem et al. 2011), and the incidence of strong rotation 
signatures in quiescent galaxies does not change from z~ 1 
to the present day (at ~ 60%, van der Wei & van der Marel 
2008). The picture that emerges here is that the vast ma- 
jority of quiescent systems have undergone some event that 
both steepens their light profile (either by creation of a dis- 
tinct bulge component, or simply by steepening the light pro- 
file) but manages to retain a significant fraction of the sys- 
tem’s original angular momentum in preserving a disk com- 
ponent. In the context of galaxy merging, such systems are a 
relatively natural outcome of merging between disk galaxies 
with even modest gas fractions, where higher mass ratio mi- 
nor mergers (or mergers between more gas-rich systems) lead 
to are progressively more disk-dominated remnants (see, e.g., 
Naab et al. 2006, Hopkins et al. 2009, Hoffman et al. 2010). 

Secondly, a small fraction of quiescent systems appear 
to have little in the way of a bulge component (Bell 
2008; Stockton et al. 2004; McGrath et al. 2008; Bundy et al. 
2010); the examples in McGrath et al. (2008) have been very 
carefully documented. At least a few of the examples in Fig. 9 
appear to have genuinely low-concentration light profiles, and 
at least one is an inclined disk. Given that such systems in the 
local universe appear to all be satellites in high mass (group 
or cluster) halos (Bell 2008), it will be interesting to explore 
the environments of such galaxies, as a function of redshift, 
and help to elucidate the extent to which disk-only quiescent 
galaxies are the products of stripping of their cold gas content 
by hot gas in a deep potential well (in this context, it is worth 
noting that the z ~ 1.5 systems in McGrath et al. 2008 were 
chosen to be in fields near radio- loud z~ 1.5 QSOs, and may 
in fact reside in overdensities). 

5.3. On the nature of high S6rsic index star- forming galaxies 

Figs. 6 and 8 show a population of star-forming galax- 
ies with high S6rsic index. We show examples of such sys- 
tems in Fig. 10 in 5 different redshift bins. One can see that 
these high $6:sic index, star-forming systems fall broadly into 
two types of systems; . At lower redshifts, there are some 
clear examples of disks with very prominent bulges; such 
galaxies are the star-forming counterparts of the quiescent 
disks discussed above, and appear to be systems that have 
either retained, or ‘re-grown’, a substantial disk during/after 
a bulge formation event (Baugh et al. 1996; Kannappan et al. 
2009; Benson & Devereux 2010). At all redshifts, there are 
what appear to be genuinely spheroidal galaxies, but in many 
cases with significant asymmetries, and in some cases evi- 
dence of interactions (e.g., nearby companions, tidal tails). 
It will be interesting to examine such systems in the future 
to understand the degree of overlap between these systems 
and ‘blue spheroids’ (Menanteauet al. 2001; HauBler 2007; 
Kannappan etal. 2009; Gyory&Bell 2010), post-starburst 
galaxies (Vergani et al. 2010), or ongoing ‘bulge-building’ 
starbursts (Wuyts et al. 201 lb) and to ask if the properties of 
these galaxies are more consistent with an interpretation as a 
quiescent galaxy in formation, or a galaxy in the early stages 
of re-growing a disk (Kannappan et al. 2009). 

5.4. Future steps 

There are a number of possible future directions for pursu- 
ing this line of study further. Sharpening the analysis would 
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Fig. 10. — F160W Postage stamps of star-forming galaxies with n > 3 and 
stellar surface densities between 1 0 9 and 1O 1O M0 kpc -2 . The postage stamps, 
within each class, are ordered by stellar mass (ordered left to right). At all 
redshifts, all postage stamps are 40 physical kpc on a side, and are scaled 
to a constant ’stellar mass density’ (total intensity is scaled to total stellar 
mass, meaning that if the stellar M/L is constant over the face of the galaxy 
this postage stamp should reflect the stellar mass density) and are displayed 
using asinh scaling (linear at low intensity, and logarithmic at higher surface 
brightness; Lupton et al. 1999). 

be possible with the addition of reliable bulge/disk decompo- 
sitions: the hypothesis that the existence of a prominent bulge 
is necessary but not sufficient to shut off star formation can be 
better tested, and additional parameters can be explored for 
possible relevance (e.g., bulge/total ratio, bulge Sdrsic index, 
or bulge mass; see, e.g., Drory & Fisher 2007 for an explo- 
ration of some of these issues with a low-redshift sample). 
More detailed exploration of the outliers at all redshifts (the 
quiescent lower Sersic index systems, or star-forming systems 
with high S6rsic index) may help to illuminate the processes 
that remove (or keep out) cold gas from galaxies. Finally, 
increasing the number statistics at z < 1 (already underway 
with, e.g., the GEMS, AEGIS or COSMOS datasets; Rix et al. 
2004, Davis et al. 2007, Scoville et al. 2007) and at z> 1 with 
the full five field coverage of CANDELS (Grogin et al. 2011; 
Koekemoer et al. 201 1) may prove useful. 

6. CONCLUSIONS 

Motivated by a desire to empirically explore the evolution- 
ary factors that lead to a lack of star formation in galaxies, 
we have explored the structures and star formation activity 
of the galaxy population using the HST/WFC3 F160W imag- 
ing data from CANDELS in the UDS field. We used pub- 
lic photometry and photometric redshifts from Williams et al. 
(2009), and determined rest-frame absolute magnitudes and 
stellar masses using our own stellar population model fits. We 
supplement this with public 24/im data from SpUDS, and sep- 
arate galaxies into quiescent and star- forming using a combi- 
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nation of optical-near-infrared colors and 24/xm information. 
For structural information, we use parametric fits of a single 
Sersic profile to the WFC3 F160W imaging data. We then 
combine these data for z > 0.6 galaxies with a sample with 
similar parameters from the SDSS, to create a z< 0.05 com- 
parison sample. We then proceed to explore the evolution of 
the A/. > 3 x 10 10 Af© galaxy population from z= 2.2 to the 
present day. 

We first visualize the evolution of the galaxy population 
over the last 10 Gyr by normalizing the sample to a ‘fixed’ 
comoving volume of 10 5 Mpc 3 . In agreement with a large 
number of other works, we find that the number density of 
massive galaxies increases approximately five-fold from z~ 2 
to the present day, and that die number density of quiescent 
galaxies increases yet more rapidly. Furthermore, examining 
the properties of the quiescent galaxy population, we find that 
the vast majority of those quiescent galaxies have high S6r- 
sic indices, a sign that they have a prominent bulge compo- 
nent. The growth of the quiescent galaxy population appears 
to be intimately linked with the growth of the high S6rsic in- 
dex population. 

Given the rapid evolution of the quiescent galaxy popula- 
tion, we proceed to explore the strength of the correlation be- 
tween quiescence and four galaxy/structural parameters: stel- 
lar mass, ‘velocity dispersion* M/R oc cr' 2 , stellar surface den- 
sity and SSrsic index. At all redshifts z< 2.2, stellar mass cor- 
relates poorly with quiescence. Bearing in mind the < 0.2 dex 
scatter between stellar mass and dark matter halo mass at 
z ~ 0, this argues against the notion that halo mass alone is 
the main determinant of quiescence. 

At z< 0.05, we find that Sersic index correlates much more 
strongly with quiescence than either velocity dispersion or 
surface density. At z > 0.6, we find that velocity dispersion, 
surface density and Sersic index correlate well with quies- 
cence, where the correlation of S6rsic index with quiescence 
is marginally stronger. 

All correlations have substantial scatter, however. Many 
quiescent systems have prominent disks, although the vast 
majority of quiescent galaxies with disks have prominent 
bulges. A very small fraction of quiescent galaxies appear to 
be disk-dominated; in the local Universe they are all satellite 
galaxies in galaxy groups/clusters, at higher redshifts we did 
not explore environmental variables for lack of dynamic range 
in galaxy environments. Star-forming systems with high n , M, 


E and o' are not particularly uncommon; at high redshifts they 
appear to be genuinely compact with high n, and often show 
asymmetries or signatures of tidal interactions (one may wish 
to associate these with the possible remnants of gas-rich dis- 
sipational galaxy interactions/mergers), and at lower redshifts 
there is a mix of similar systems and composite bulge-disk 
star-forming systems (with large bulges). 

At z< 0.05, Bell (2008) concluded that a prominent bulge 
(and by association, a supermassive black hole) was a neces- 
sary but not sufficient condition for a galaxy to turn off its own 
star formation on galaxy-wide scales (all quiescent galaxies 
in the centers of their own halos had prominent bulges, but 
not all galaxies with bulges lack star formation). This obser- 
vational association is qualitatively consistent with the AGN 
feedback paradigm (no supermassive black hole, no ability to 
shut off star formation). While there is clearly scope for fur- 
ther investigation of the drivers of quiescence, the evidence 
assembled here appears to be consistent with this proposition 
to z< 2.2, a time interval of more than lOGyr. 
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APPENDIX 

APPENDIX A: NORMALIZING THE SAMPLE TO A CONSTANT COMOVING VOLUME 

In Fig. 3, we Monte Carlo subsampled the galaxy population of the UDS and the SDSS to an equivalent comoving volume of 
lC^Mpc 3 . There are two ways to do this: simply rescaling the sample by 10 5 Mpc 3 / V(z), where V(z) is the comoving volume 
of that redshift bin (leaving one susceptible to the first-order effects of sample variance), or by rescaling the sample in a given 
redshift bin to have the number of galaxies expected at that redshift above that mass limit using mass functions derived from 
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FIG. 11. — The number of galaxies with stellar masses in excess of 3 x 1O 1O M 0 in a volume of K^Mpc 3 measured from the UDS (diamonds, with the dominant 
sample variance uncertainties as estimated using the method of Somerville et al. 2004; see also Moster et al. 2011), and the number of galaxies with masses in 
excess of 3 x 1O 1O M 0 (crosses) inferred from the stellar mass functions of Ilbert et al. (2010, z > 0.2) and Bell et al. (2003, local values). The number of galaxies 
with M* > 3 x 10 10 Mq in the UDS is consistent to within the significant sample variance and systematic stellar mass uncertainties with those of larger fields. We 
choose ir. this paper to adjust the number of galaxies for display in Fig. 3 to fit the straight line in this figure, to take out the first-order effects of sample variance 
on the results (while noting that doing so makes no important difference to our conclusions; the evolution of the population is not subtle). 


much larger surveys (canceling out number density variation from sample variance, but leaving behind any systematic variation 
in galaxy properties that are a function of the variation in the average environment in this light cone as a function of redshift). 

We compare the two approaches in Fig. 11. Diamonds with error bars denote the number of galaxies with stellar masses in 
excess of 3 x 1O 1O M 0 in the UDS, with error bars denoting the expected degree of sample variance in that bin following the 
method of Somerville et al. (2004). Crosses at z> 0.2 show the number of galaxies with stellar mass in excess of 3 x 1O 1O M 0 
expected from the stellar mass functions of Ilbert et al (2010). The cross at low redshift shows the number of galaxies with 
M* > 3 x 1O 1O M 0 from the mass function of Bell et al. (2003). The number of galaxies observed in the UDS is broadly consistent 
with, or perhaps somewhat larger than, the number of galaxies expected from larger surveys, given the substantial sample variance 
uncertainties. There are systematic uncertainties also on the crosses; choosing different star formation histories for constructing 
stellar masses, stellar population models, etc., can give more than a factor of two variation in stellar masses that translates 
into around 50% number density uncertainty. For Fig. 3, we chose to rescale the number of galaxies to the smoothly- varying 
number of galaxies given by the line shown in Fig. 11 (approximately corresponding to scaling the number of galaxies to larger 
cosmological surveys). 

None of the results shown in Fig. 3 depend in any important way on the choice of this scaling method; the evolution of the 
galaxy population from z ~ 2 to the present day is not subtle and is robust to even the significant systematic uncertainties inherent 
to mass function analyses. While this scaling in the number of galaxies largely counteracts the worst effects of sample variance, 
if the properties of galaxies depend strongly on environment there will be a second-order difference between the properties of 
galaxies sampled from a globally-underdense vs. a globally-overdense volume. In the absence of enough environmental dynamic 
range to robustly measure its effect within CANDELS at the current time (see Papovich et al. 2011 for a discussion of SFR 
and galaxy size differences as a function of environment of this using a known z- 1.62 galaxy cluster in the CANDELS UDS 
coverage), we made no attempt to correct for this second-order effect in this paper. 




